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I. INTRODUCTION 

The emergence of classical dynamical chaos from 
more profound quantum mechanics is one of the most 
intriguing problems in physics h, _2., _4] . Dynamical 
chaos in classical mechanics is a special kind of ran- 
dom motion in dynamical systems without any noise 
and random parameters that is characterized by sen- 
sitive dependence on initial conditions in a bounded 
phase space. The current consensus is that isolated 
bounded quantum systems do not show sensitive de- 
pendence on initial conditions in the same way as clas- 
sical systems because their evolution is unitary. So the 
question is: What is the fundamental mechanism of 
arising classical chaos from quantum mechanics? 

The purpose of this lecture is to show how classi- 
cal Hamiltonian chaos with sensitive dependence on 
initial conditions, positive Lyapunov exponents, frac- 
tal properties of underlying phase space, and anoma- 
lous statistical characteristics may arise from quan- 
tum electrodynamics of a single atom strongly inter- 
acting with a quantized mode of the electromagnetic 
radiation field in a high-quality cavity. The study of 
the fundamental atom-photon interaction constitutes 
the rapidly growing field of cavity quantum electrody- 
namics (for a review see 

I1I10I3)- The experimental 
state of the art has reached in this field the stage where 
the transition from classical to quantum dynamical 
regimes can now be probed directly. Atoms and pho- 
tons, confined in a high-quality cavity, are ideal ob- 
jects to study quantum-classical correspondence and 
quantum chaos. 

Real quantum systems are not isolated. They in- 
teract with their environment and in attempts to 
measure their states with classical measuring devices 
which should, by virtue of their purpose, be in unsta- 
ble states. So we are dealing with quantum-classical 
hybrids. It is not a simple question which degrees 
of freedom should be considered in a given physical 



situation as quantum ones and which ones as clas- 
sical. Let us consider an excited atom in a single- 
mode high-quality cavity whose frequency is close to 
the frequency of one of the atomic electrodipole tran- 
sitions. The atom emits a photon in the cavity mode 
and goes to a lower-lying state. Under the conditions 
of the strong atom-field coupling, the photon may be 
reemitted and reabsorbed by the atom many times. 
Sooner or later, the atom will emit a photon in one of 
the other modes of the electromagnetic field that are 
not sustained by the cavity (for example, transverse 
modes in an open Fabry-Perot cavity) and the photon 
will be lost in a surrounding, the process known as 
spontaneous emission. In real experiments, the cavity 
relaxation due to the losses in the cavity walls cannot 
be neglected. It is usually modeled by a weak coupling 
of the selected cavity mode with a bath of harmonic 
oscillators (spanning a wide frequency range) in ther- 
mal equilibrium at a given temperature. Sooner or 
later, the atom will relax to a lower-lying state that 
is not resonant with the cavity mode. It seems to 
be reasonable to treat two near-resonant atomic lev- 
els, strongly coupled to a cavity mode, as a quantum 
dynamical system weakly coupled to the other elec- 
tromagnetic modes, the other atomic levels and the 
cavity walls. An infinite number of the respective de- 
grees of freedom, which are coupled to the quantum 
dynamical system's degrees of freedom but are not af- 
fected by them, forms an inexhaustible external reser- 
voir (that may be treated classically) usually called 
"the environment". It is a kind of the external cou- 
pling which is inevitably present in reality. 

There is another kind of coupling called the dynam- 
ical coupling. For example, if the average number of 
photons in the cavity mode is sufficiently large, one 
may treat the field as a classical object dynamically 
coupled to quantized atoms (the semiclassical approx- 
imation) and take into account the feedback effect of 
the atoms on the radiation field. The semiclassical 
approximation breaks down when we deal with sin- 
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gle cold atoms and photons in a cavity. In this sit- 
uation we should take into account the translational 
(external) atomic degree of freedom. When a cold 
atom emits and absorbs photons its momentum and 
position may vary significally due to the recoil effect. 
The external atomic degree of freedom, which may 
be treated classically if the values of the atomic mo- 
mentum are greater than the photon momentum, is 
dynamically coupled to the internal atomic and field 
degrees of freedom which are treated quantum me- 
chanically. We show in this lecture that, even in the 
absence of any other interaction with the environment, 
the interaction of the quantum atom-field system with 
the external atomic degree of freedom provides, un- 
der appropriate conditions, the emergence of classical 
chaos from quantum electrodynamics. 



II. EARLY STUDIES OF CLASSICAL CHAOS 
IN THE ATOM-FIELD INTERACTION 

The discovery that a single-mode laser, a symbol 
of coherence and stability, may exhibit deterministic 
chaos is especially important not only because lasers 
are one of the main instruments in physics but lasers, 
as well, provide almost ideal systems to test general 
ideas in quantum mechanics and statistical physics. 
From the standpoint of nonlinear dynamics, laser is 
an open dissipative system that transforms an exter- 
nal excitation into a coherent output in the presence 
of loss. Some manifestations of strange attractors and 
dissipative chaos have been observed with different 

types of lasers (for a review, see inaiiii). A col- 

lection of identical two- level atoms, interacting with 
a single-mode electromagnetic field, provides the sim- 
plest model for laser dynamics. Because of a large 
number of atoms (and photons), laser dynamics can be 
adequately described in the so-called semiclassical ap- 
proximation, where one treats atoms as two-level ob- 
jects, which may be described by the Bloch equations 
for two components of the collective electrodipole po- 
larisation and collective atomic population inversion, 
interacting with a field mode to be governed by the 
classical Maxwell equations for the field strengths with 
the right hands depending on the atomic polarisation. 
With a single-mode homogeneously broadened laser, 
operating in resonance with the gain center, these five 
equations can be reduced to three real-valued equa- 
tions for slowly varying amplitudes which have been 
shown to be equivalent to the well-known Lorentz 
model for fluid convection [l^ . 

Practically in the same time, ideas of dynamical 
chaos have been explored with fundamental models 
of the matter-radiation interaction, comprising of a 
collection of two-level atoms interacting with their 
own radiation field in a perfect single-mode cavity 
without any loss and external excitation. It has 
been shown theoretically and numerically 1^ that 
the semiclassical Maxwell-Bloch equations, following 



from the Dicke Hamiltonian |14| . may demonstrate 
Hamiltonian semiclassical chaos if one goes beyond 
the so-called rotating-wave approximation, i. e. if 
one takes into account energy non-conserving terms 
in the Hamiltonian (for details, see the beginning of 
the next section). This mechanism of arising Hamil- 
tonian semiclassical chaos is rather weak under real- 
istic assumptions. On the other hand, rotating-wave 
approximation suppresses chaos with atoms at rest 
due to existence of an additional integral of motion, 
the conserving interaction ener gy. It has been shown 
in a series of our papers [isl Il6l 11^ EM that 
Hamiltonian chaos may arise within the rotating-wave 
approximation under conditions of a modulation of 
the atom-field interaction. In a natural way it oc- 
curs when atoms move through a cavity in a direc- 
tion along which the cavity sustains a standing-wave 
field that is periodic in space 13- When mov- 
ing with a constant velocity, atoms "see" the field 
whose strength changes in time periodically. It breaks 
down the interaction-energy integral and may cause 
Hamiltonian semiclassical chaos of a homoclinic type 
0,0,0]. There are another ways of modulations to 
be considered in [tR |25 | . The modulation of the de- 
tuning between the atomic transition frequency and 
the frequency of the field mode has been shown to 
produce parametric instability and Hamiltonian semi- 
classical chaos JJ7] . Structural Hamiltonian chaos 20| 
may arise as a result of a harmonic modulation of cav- 
ity length which causes the respective oscillations of 
the nodes of a standing-wave. 

Generally speaking, the model of the atom-field in- 
teraction should involve not only the internal atomic 
and field degrees of freedom but also the center-of- 
mass motion of the atom. When emitting and ab- 
sorbing photons, atoms not only change their internal 
states but their position and momentum are changed 
as well due the photon recoil effect. This effect may 
be neglected if one deals with Rydberg atoms inter- 
acting with a microwave field in a cavity (as it has 
been done in [H, [S El El H^) or with thermal 
usual atoms interacting with a visible light. It has 
been theoretically and numerically shown in [2ll |2^ 
that Hamiltonian semiclassical chaos may arise with 
cold atoms with recoil in a standing- wave microcavity. 

In trying to describe adequately dynamics of single 
atoms and photons in a high-quality cavity, one should 
go beyond the semiclassical approximation and treats 
the atom-photon interaction on a quantum ground. 
The fully quantum model of the interaction between a 
single two-level atom (without a recoil) with a single- 
mode quantized field in an ideal cavity is known as 
the Jaynes-Cummings model j^^- It describes the 
atom-field system as a quantum-electrodynamical ob- 
ject whose evolution in time is (quasi)periodic. In 
Schrodinger picture, it may be described by an in- 
finite set of linear ordinary differential equations for 
the probability amplitudes to find the atom in the 
ground/ excited state and the field in the state with 
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n photons, where n runs from zero to infinity. When 
adopting the semiclassical approximation, we reduce, 
by hook or by crook, this infinite set to a smaU num- 
ber of equations for atomic and field variables. In 
fact, we decouple the quantum atom-field system into 
an atom and field parts which may exchange excita- 
tions with each other. In doing so, we get automati- 
cally products of the atom and field variables in the 
equations of motion. The respective Maxwell-Bloch 
equations may, under appropriate conditions, produce 
chaos in the classical sense of sensitive dependence on 
initial conditions in the reduced classical phase space 
spanned by the atomic and field expectation values. It 
resembles the procedure of deriving the famous Lorenz 
equations from an infinite hierarchy of mode equations 
in fluid convection by reducing it to three main modes 
only. 

The reduction of an infinite set of linear equations, 
that are not chaotic in the classical sense, to a finite 
set of nonlinear equations, that may be chaotic, makes 
an impression that it is only a mathematical trick. It 
seems that it is not only a useful trick enabling us to 
handle with the equations of motion but the reduc- 
tion is a model of processes that may occur in nature. 
Coherence loss caused by inevitable interaction with 
environment (decoherence) breaks quantum unitarity 
suppressing some quantum properties of motion and 
revealing its classical properties. Such a situation can 
be modeled with the use of a quantum-classical hy- 
brid, a system with quantum and classical degrees of 
freedom dynamically coupled to each other. Which 
part of the whole system under consideration is quan- 
tized and which one is treated classically depends on 
the physical situation. If we have a large number of 
atoms (or photons) we may adopt with a good ac- 
curacy the semiclassical approximation treating the 
cavity field as a classical one. When considering a 
single atom with a recoil in a high-quality microcav- 
ity it is necessary to quantize the internal electronic 
atomic and field degrees of freedom but the transla- 
tional atomic motion may be treated classically if the 
values of the atomic momentum are greater than the 
photon momentum. 



III. ATOM WITHOUT RECOIL 
INTERACTING WITH A QUANTIZED 
RADIATION FIELD IN AN IDEAL CAVITY 

The interaction between matter and radiation, com- 
monly evidenced by spontaneous emission and in- 
duced emission and absorption of photons by atoms, is 
one of the most fundamental dynamical interactions in 
nature. In free space, an excited atomic state decays 
irreversibly because an infinity of vacuum electromag- 
netic states is available to the radiated photons which 
disappears in the free-space mode continuum. The sit- 
uation may be changes if the mode structure and the 
density of modes are modified by placing the radiating 



atom into a high-quality cavity. A cavity-induced en- 
hancement and inhibition of spontaneous emission has 
been observed in a number of experiments in different 
ranges of the electromagnetic fieldjfrom microwaves 
to visible light (for a review see |a.l24l|). 

In a near-resonant ideal cavity with a single-mode 
frequency closed to the atomic transitional frequency, 
the radiation emitted by the atom is reflected at the 
cavity walls and is reabsorbed by the atom many times 
before it dissipates. In particularly, spontaneous emis- 
sion becomes reversible. In cavity quantum electrody- 
namics near-resonant interaction of a single atom with 
a single-mode cavity field is commonly modeled by the 
Hamiltonian 

H = Ha + Hf + Hint, (1) 

where 

Ha = ^huJa^z (2) 

represents the free motionless two-level atom with fkUa 
being the energy separation between the electronic 
states, the ground |1) and excited |2) ones. The Pauli 
spin operator acts on the states as 

a,|l> = -|l), a,|2)-|2). (3) 

The Hamiltonian 

Hf ^ n^f{a^a + 1) (4) 

represents the mode of the isolated field. The photon 
annihilation a and creation a) operators with the com- 
mutation rule [a, a)] = 1 act on the photon number 
states as follows: 

a|n) = Vn|n- 1) , a) \n) = %/n + 1 |n + 1) . (5) 

The electric dipole interaction between the atom and 
the field mode is described by the operator 

Hint = hVtQ{x){a + a)){cr+ + ct-) ~ 

Ml^){x){ad-j^ + a) a_)^ (6) 

where ^q{x) is the atom- field coupling constant at 
the atomic position x in the one-dimensional cavity, 
which is known under the names "vacuum or single- 
photon Rabi frequency" . The Pauli spin operators 
CT± with the commutation rules [^±,(72] — t2<5'±, 
[(7+, a-] = az describe transitions between the ground 
and excited electronic states 

a+|l> = |2), a_|2> = |l>, a+|2)=a_|l>-0. 

(7) 

The first form of the Hamiltonian ^ follows from 
the quantization of the classical electric dipole atom- 
field interaction by replacing the atomic dipole mo- 
ment by the operator = cr^ + (7_ and the classical 
electric field by the operator a -\- cl^ . The products 
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a^(T_ and aa^ describe the transition of the atom from 
the ground (excited) to the excited (ground) state 
and simultaneous annihilating (creating) a photon in 
the field mode, respectively. They are usual energy- 
conserving processes. The product describes the 
transition of the atom from the ground to the excited 
state and simultaneous creating a photon in the mode. 
The reverse process of simultaneous atom and field de- 
excitation is described by the term aa- . The last two 
processes are not energy-conserving and may be ne- 
glected while we deal with a single atom. Omitting 
the energy-nonconserving terms defines the Jaynes- 
Cummings model j23| which was proposed originally 
in 1963 as a purely theoretical tool for studying fun- 
damentals of the interaction between matter and radi- 
ation. It seemed to be far from reality because neither 
the atom nor the cavity mode can hardly be perfectly 
isolated from environment with those-days technique. 
Exciting progress in experimental techniques has dras- 
tically changed the situation. Experiments, using very 
high-quality microwave cavities with Q ~ 10^° and op- 
tical microcavities with Q ^ 10®, have now achieved 
the exceptional circumstance of strong coupling be- 
tween atoms and a cavity field with the strength of the 
coupling exceeding the atomic and cavity decays that 
provides manipulations with single atoms and photons 
and experimental provingof some predictions of the 
Jaynes-Cummings model The art of exper- 

imentalists and theoretical efforts have opened excit- 
ing perspectives for realizing gedanken experiments on 
fundamentals of quantum mechanics and implement- 
ing quantum communications and computing. As we 
will show in this lecture, an atom, interacting with a 
radiation field in a high-Q cavity, provides an exam- 
ple of strongly- coupled microscopic nonlinear system 
whose dynamics may be very complicated and even 
chaotic. 

The Jaynes-Cummings Hamiltonian has the form 

Hjc = ^hwacr^ -\- hbjf{a^a + i) + 

h^of{x){a)a-+ad+), (8) 

where is an amplitude value of the vacuum Rabi 
frequency and f{x) is a shape function of the cav- 
ity mode. Let us expand the atom-field state vector 
over the electronic atomic states and Fock (or photon- 
number) field states \n) 



results in an infinite-dimensional set of coupled linear 
equations for the complex-valued probability ampli- 
tudes 



-i [Acan + \/n+ I f{x)b. 



n+l \ , 



bn+i = i [Af,6*_,_;^ + Vn+ l/(a;)a*] 



(11) 



where dot denotes differentiation with respect to di- 
mensionless time r = ^^t, — {nujf + uia/'2)/i}Q, 
and Afc = [{n + l)w/ — Wa/2]/ilo- It is convenient to 
define new real- valued combinations of the probability 
amplitudes 



2Re(a„6;+i), 



Vn 



-2Im(a„6,*+J, 



for which we get from Hll|) the quantum Bloch-like 
equations of motion 



iln = 6Vn, 

Vn = ~SUn - 2\/rt + 1 !{x)Zn, 



(13) 



i„ = 2\/n + l ,f{x)vn, n = 0, 1, 2, 



where 6 — {uja — f)/^o is the dimensionless detuning 
between the atomic transition, Ua, and the field- mode, 
LOf, frequencies. For each specified photon number n, 
the quantity Rf^ = u\ -\- + is conserved. Two 
global integrals of motion 



Wjc = ^ Vn + 1 f{x)un - T^^^n, 



^i?„ = l (14) 



n=0 

reflect conservation of the total energy and of the to- 
tal probability, respectively. Eas. ((T^ can be easily 
solved for each n in terms of trigonometric functions 
(we assume that f{x) does not depend on time). The 
general solution of the Schrodinger equation is the sum 
of the solutions of an infinite set of the independent 
Bloch-like equations. The exact general solution for 
one of the measured quantities, the atomic population 
inversion 



oo 

z{t) = ^z„(t), 

n=0 



(15) 



>„(i)|2,n) +6„(t)|l,n) 



(9) 



is the following: 



where a„(6„) are probability amplitudes to find the 
atom in the excited (ground) state and n photons in 
the field mode, respectively. Substitution of Eqs.® 
and into the Schrodinger equation 



ih 



d\^{t)) 
dt 



Hjc !*(<)) 



(10) 



z(r) 



2^W„(0) -2 (1 -COSfi, 



, ,2v^m/(x) . ^ 

(16) 
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where 



(a) 



v/<52+4(n + l)/2(x) 



(17) 



is the n-photon Raby frequency, m„(0), w„(0) and 
z„(0) are the respective initial values which are de- 
termined by the initial atomic and field states. As it 
is seen from (I16|) . the solution with initially sharply 
defined atomic and field energy states (|1) or |2) and 
\n)) describes a periodic exchange of one quantum of 
energy between the atom and the field. 

The Fock state \n) with a specified number of pho- 
tons in the mode is an exotic field state. In general, a 
pure quantized field state is an infinite superposition 
of the photon-number states 



I/) = ^Cn \n) 



(18) 



n=0 



where p„ 



is the probability for observing n 



photons. The most classical of single-mode quantum 
states is a coherent state 



(19) 



n=0 



whose photon probabilities follow a Poisson distribu- 
tion with the mean number of photons < n >= jap 
and the root mean square spread An = Evolution 
of the population inversion of an atom, z, in the field, 
initially prepared to be in a coherent state, contains 
all the Rabi frequencies 0„. The incommensurability 
of Rabi frequencies for different n inevitably washes 
out the periodicity of population transfer resulting in 
a collapse of the excited-state occupation probability. 
Since the frequencies r2„ form a discrete set, the signal 
reappears after a time (see Fig. QJi computed with 
f{x) — I, S ~ and with the initial conditions (j^ ). 
The collapse-revival phenomenon, a direct demonstra- 
tion of quantum nature of the radiation field, has been 
demonstrated experimentally 0, . 

Let us consider as well another pure field state 



(20) 



n=0 



whose photon probabilities follow a Bose-Einstein dis- 
tribution 



Pnix) = 



< n >' 



(1+ < n >)"+!■ 



(21) 



Collapses and revivals also occur with this field state 
(see Fig. ^ computed with f{x) = 1, S = and with 
the initial conditions H36(l 'l. but the collapse time is 
much shorter as compared with a coherent state be- 
cause the Bose-Einstein spread in the photon number 
is far larger than for a coherent field with the same 
mean number of photons. 

Fig. n demonstrates that the simple Jaynes- 
Cummings model with f{x) = 1 may produce very 
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FIG. 1: Atomic population inversion z{t) for the resonant 
Jaynes-Cummings model with an initially excited motion- 
less atom and (a) coherent and (b) Bose-Einstein initial 
states of the cavity quantized field. 



complicated time evolution. However, the respective 
oscillations, of course, are not chaotic in the classi- 
cal sense of exponential sensitivity to small changes 
in initial conditions. It is clear from the very struc- 
ture of the equations of motion H13|) which are infinite- 
dimensional but linear. In classical mechanics the 
phase space is continuous and states of a classical sys- 
tem can be arbitrarily close to each other providing 
a possibility of exponential diverging of initially close 
trajectories in a bounded region of the phase space. 
In quantum mechanics there is no notion of a trajec- 
tory, and the quantum phase space is not continuous 
due to the Heisenberg uncertainty principle. Evolu- 
tion of an isolated quantum system is unitary and no 
dynamical chaos in the classical sense is possible with 
such systems. The notion "quantum chaos" just refers 
to the behavior of quantum systems whose classical 
counterparts behave chaotically 0, Real quantum 
systems are not isolated, they interact with their envi- 
ronment and, under attempts to measure their states, 
with classical devices which, by virtue of their pur- 
pose, should be in unstable states. Loosing of co- 
herence, due to inevitable interaction with environ- 
ment (decoherence), breaks down quantum unitarity 
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suppressing quantum properties of motion and mani- 
festing classical ones. A quantum rotor with periodic 
kicks provides an illustrative example of suppressing 
quantum dynamical localization when one takes into 
account a thermal "bath" in the respective equations 
of motion In the next section we will elaborate 
these ideas with an extended version of the strongly 
coupled atom-field system considered above. 

IV. A QUANTUM-CLASSICAL ATOM-FIELD 
HYBRID 

In the process of emitting and absorbing photons 
of the cavity-field mode, the atom not only changes 
its internal electronic states but its external trans- 
lational state is changed as well due to the photon 
recoil effect. In this section we consider a single two- 
level atom with mass nia moving in an ideal cavity 
which sustains a single standing-wave mode along the 
axis X with the wave vector k f and the shape function 
f{x) = — cos{kfx). The respective Hamiltonian is the 
following extension of the Jaynes-Cummings Hamilto- 
nian ( IHl) 

H = p + —huJaO'z + hcufa'a— 

fi,ilo (a^fT- + ao--|_) coskfx, (22) 

where the momentum p and position x operators sat- 
isfy the standard commutation relation [x,p\ — ih. 
Operators, belonging to different degrees of freedom, 
commute with each other at the same time moment. 

We have now three degrees of freedom, the internal 
atomic and the field ones, and the external (or trans- 
lational) atomic degree of freedom. The first two de- 
grees of freedom are treated as fully quantum ones in 
the Schrodinger picture. In fact, there are an infinite 
number of quantum degrees of freedom (see Eq. 
which are entangled. The external degree of freedom 
will be treated as the classical one that may be justi- 
fied by large values of the atomic momentum as com- 
pared with the photon momentum hkf . The Hamilton 
equations of motion for the classical external degree 
of freedom is easily found from the Hamiltonian (|22ll 

d<x>_d<H> d<p>_ d<H> 

dt d <p> ' dt d < X > ' 

(23) 

where < • • • > denotes an expectation value of the 
corresponding operator over a quantum state j^*) of 
the electronic-field Hamiltonian. Using the normal- 
izations X = kf < X >, p =< p > /hkf, and r — Slot, 
we get from 

X — ap, 

I ,.t^ I \ . (24) 

p = — \^'(''") a^(T_ + a(j_|_ ^{t)) sinx, 

where a — kkJ/niai^Q is the normalized recoil fre- 
quency which characterizes the average change in ki- 
netic energy of the atom, h^k'l/2ma, in the process of 



emission and absorption of a photon. After computing 
the expectation value in H24|) with the state vector @ , 
we obtain our basic Hamilton-Schrodinger equations 



X 


= ap, 




oo 


V 


= — V^i + 1 it„ sin X, 




n=0 


Un 


= SVn, 


Vn 


— —SUn + 2\/?l + 1 2„ cos X, 




= —2y/n + Ivncosx, n 



(25) 

= 0, 1, 2, ... 



which describe the atom-field quantum-classical hy- 
brid with the quantum degrees of freedom dynami- 
cally coupled to the classical degree of freedom. This 
infinite set of nonlinear ordinary differential equations 
possesses an infinite number of the integrals of motion, 
the total energy integral 

oo S ^ 

W = - - ^ \/n + 1 -u„ COS x - - ^ 2„, (26) 

Tl = 11=0 

the Bloch-like integral for each n 

Rl = < + < + zl, (27) 
and the global integral 

oo 

Y^Rn^l- (28) 

n=0 

The infinite-dimensional nonlinear dynamical sys- 
tem (|25|l is a quantum generalization of the five- 
dimensional semiclassical set of the equations of mo- 
tion for the same problem that was derived in Refs 
l2U|^ with a classical field. The later one was shown 
|2ll I22 . 26, 27j to be chaotic with positive values of 
the maximal Lyapunov exponent in some ranges of 
the system's control parameters. It is rather difficult 
to compute this quantitative indicator of dynamical 
chaos with an infinite dimensional set of ODE's. How- 
ever, we have found another signatures of Hamiltonian 
chaos with Eqs. H25|l . the main of which are fractals 
which will be demonstrated in the next section. 

Let us compare the Rabi-oscillation signals with the 
atom-field quantum-classical hybrid H25I) and with the 
fully quantum Jaynes-Cummings model without recoil 
H13I) . In Fig. [2Ji and b the population inversion z (see 
Eq. H15|) ') for an initially excited atom interacting in 
resonance with the field that is initially in a coherent 
state (|19|l and in a Bose-Einstein state (|20|l . respec- 
tively, is shown. This figure should be compared with 
Fig. n] where the Rabi oscillations with the Jaynes- 
Cummings model have been computed for the same 
initial conditions and the field states with < n >= 10 
but with a motionless atom. More pronounced col- 
lapses and revivals occur with the quantum-classical 
hybrid as compared with the fully quantum model 
without recoil. In order to understand peculiarities of 
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FIG. 2: The same as in FIG. Q but with the resonant 
quantum- classical hybrid. 



the Rabi oscillations in Fig. El look at Fig.[3i demon- 
strating the scheme of gedanken experiments. An 
atom starts at x = and moves to the right with 
the initial momentum pq = 50. It reaches the first 
node of the standing wave, where the coupling coef- 
ficient with the field mode is zero at time moment 
Ti — 7r/2ckpo — 31.4 in dimensionless units (we choose 
a — 0.001 in our computer simulations). For both 
the field states, one can see the respective slowing 
down of the Rabi oscillations at the moments of time 
Ts — (l + 2s)T:/2apQ, when the atom transverses the s- 
th node. The pronounced peaks of revivals of the Rabi 
oscillations occur at the moments = rir/apo, when 
the atom transverses the r-th antinode of the standing 
wave where its coupling with the field is maximal. 

In spite of a complicated character of the Rabi os- 
cillations, shown in Fig. [3 they are regular because at 
exact resonance, 6 = 0, all the values of the real- valued 
amplitudes Un are conserved during the evolution. 
With initially excited atom we have m„(0) = w„ (0) = 
for each n. It immediately follows from Eas. 1251) that 
the atom moves with a constant velocity through the 
cavity. It means that x is a linear function of time, 
and the Hamilton-Schrodinger equations (|25|l reduce 
to the periodically modulated linear Bloch-like equa- 



tions (|13|l with periodic solutions. The exact solu- 
tion for the atomic population inversion can be easily 
found 



z{t) = ^ z„(0) cos 



n=0 



2yM + T 

apo 



■ sm apoT 



(29) 



It is a periodic function (shown in Fig. |2J) with the 
period to be equal to n/apo and the maxima at = 
rn/apQ (r = 0,1,2,...). 

Out off resonance, d ^ 0, the quantum-classical hy- 
brid may demonstrate chaos. To diagnose chaos it is 
instructive to compute the maximal Lyapunov expo- 
nent A whose values depend on initial conditions, on 
the detuning i5, the mean number of photons in the 
mode, and on the recoil frequency a. If the quantized 
field is initially prepared in the Fock state |n) with ex- 
actly n quanta in the mode, the infinite-dimensional 
set reduces to 8 equations 



X 

P 

Un-l 
Vn-l 
Zn-1 
Un 

Vn 



ap, 



/nun-i + V'^ + 1 Un) sinx, 

— 5Un-l + 2y/nZn-l cosx, 
Vn-l COS X, 

—Sun + 2\Jn + 1 z„ cosx. 



(30) 



-2\/n -I- 1 Vn 



cos X 



with initial conditions 

a;(0) = .To, _p(0)=po, 

^„-i(0) = -|fo„(0)|', z„(0) = |a„(0)|^ 

Un-M = w„(0) = Vn-M - w„(0) - 



(31) 



describing the atom-field system initially prepared in 
a superposition state 



l^'(O) >= a„(0)|2, n > +fe„(0)|l, n >, 
|a„(0)|' + |&„(0)|' = l. 



(32) 



The so-called topographic A-maps, showing by color 
modulation values of A in the ranges of values of two 
control parameters with the third one to be fixed, 
have been computed with the set H30I) . In Fig. Owe 
show the A-map with initial zero atomic population, 
a„(0) = 6„(0) = 1/V2, z{0) = z„_i(0) + z„(0) = 0, in 
the ranges of the detuning S and the recoil frequency 
a with the fixed value of the initial number of photons 
n = 10. At exact resonance, the set (|30|l is integrable 
and A = at i5 = 0. In the range a ~ 10~'*-^10~^, that 
corresponds to realistic values of the recoil frequency, 
the maximal Lyapunov exponent may be positive and 
one may expect chaotic atomic motion in the respec- 
tive ranges of a and S. Another A-map with initially 
excited atom, a„(0) = 1,6„(0) = 0, z(0) = 1, demon- 
strates in Fig. 0] the values of A in dependence on a 
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FIG. 3: The maximal Lyapunov exponent A with the 
quantum-classical hybrid in the Fock quantized field ver- 
sus the atom-field detuning S and the logarithm of the 
dimensionless recoil frequency a. 




109,0" 



FIG. 4: The double logarithmic plot of A versus a and 
initial number of photons in the Fock field n. 



and n at S — 0.5 in double logarithmic scale. The 
magnitude of the maximal Lyapunov exponent grows, 
in average, with increasing the initial number of pho- 
tons in the cavity mode. 

Computing the values of A is difficult with the cav- 
ity field initially prepared in a superposition state 
(|18|l (including coherent and Bose-Einstein states) 
that generates an infinite number of the equations of 
motion (|25|l . We have found another signatures of 
chaos. In Fig. [S] an atomic trajectory, x{t), and os- 
cillations of the atomic momentum, p{t), are shown 
to illustrate manifestations of chaos in the set (|25|l 
with the initially coherent field and the excited atom 




FIG. 5: (a) A chaotic atomic trajectory with an initially 
coherent field and (b) the respective oscillations of the 
atomic momentum. 



with < 71 10, (5 0.4, a = 0.001, and po = 25. 
A typical weakly chaotic atomic trajectory represents 
a kind of a random walking with small oscillations 
of the atom in potential wells, long ballistic flights 
of the Levy type with almost constant velocity and 
erratic turnbacks. Just the presence of Levy flights 
fluenses strongly the statistical properties of the sys- 
tem generating power-like statistical laws that have 
been extensively studied in the semiclassical approxi- 
mation in Refs.p^. l2^. Fig. |H1 presents Poincare sec- 
tions of motion in the set H25() with five different val- 
ues of the initial atomic momentum to be projected 
on the plane of the atomic external variables 
Fig. demonstrates an example of chaotic Poincare 
section {6 = 0.1) whereas Fig. [TJd presents a regular 
motion under the same conditions except for the de- 
tuning S — 0.5. 

A feasible scheme for detecting manifestation of 
chaos with hot two-level Rydberg atoms moving in 
a high-Q microwave cavity has been proposed in [2a| . 
The same idea could be realized with cold usual atoms 
in a high-Q microcavity. Consider a 2D-geometry of 
a gedanken experiment with a monokinetic atomic 
beam propagating almost perpendicularly to the cav- 
ity axis X. In a reference frame moving with a constant 
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(a) 




(b) 



FIG. 6: Projections of the Poincare sections on the plane 
of the atomic momentum p and the position x. (a) Chaos 
at J = 0.1 and (b) regular motion at 6 = 0.5. 



the formula l|33|l can be very short: with A =0.5 the 
predictability horizon Tp may be of the order of 10 in 
units reciprocal of the vacuum Rabi frequency Qq that 
corresponds to tp ~ 10~^ s with the realistic value of 
Qo ~ 10^ rad- s'^. 



"Hi 




velocity in the perpendicular direction, there remains 
only the transverse atomic motion along the axis x. 
One measures atomic population inversion after pass- 
ing the interaction zone. Before injecting atoms in the 
cavity, it is necessary to prepare all the atoms in the 
same electronic state, say, in the excited state, with 
the help of a 7r-pulse of the laser radiation. It may be 
done with only a finite accuracy, say, equal to Azjn 
for the initial population inversion Zi^. The values of 
the population inversion Zout are measured with detec- 
tors at a fixed time moment. If we would work with 
the values of the control parameters corresponding to 
the regular atom-field dynamics, we would expect to 
have a regular curve Zout--2in- In the chaotic regime, 
the atomic inversion at the output can be predicted 
(within a certain confidence interval Az) for a time 
not exceeding the so-called predictability horizon 



1 , Az 



(33) 



which depends weakly on Azjn and Az. Since the 
maximal confidence interval lies in the range |Az| ^ 1 
and A may reach the values of the order of 1.5 (see A- 
maps), the predictability horizon in accordance with 



FIG. 7: Dependence of the output values of the atomic 
population inversion Zout on its initial values Zin with S = 
0.4 (a) at r = 100 and (b) at r = 200 with the inset 
showing this dependence at the exact atom-field resonance, 
5^0. 



In the regular regime, the inevitable errors in 
preparing Azin produce the output errors Azout of the 
same order. In the chaotic regime, the initial uncer- 
tainty increases exponentially resulting in a complete 
uncertainty of the detected population inversion in a 
reasonable time. It is demonstrated in Fig. [7| where 
we plot the dependence of the values of z(t) = Zout 
at T = 100 (Fig. and r = 200 (Fig. [7|d) on the 
values of z(0) = Zjn in the chaotic regime with an ini- 
tially Fock field at S = 0.4 and A ~ 0.05. Simulation 
shows that an initial error Azin — 10^^ in preparing 
the atomic electronic state leads to complete uncer- 
tainty Azout — 2 in a rather short time. To see the 
difference, it is desirable to carry out a control exper- 
iment at the exact resonance {S ~ 0) when the atomic 
motion is fully regular with any initial values. The de- 
pendence Zout~-Zin with (S = is demonstrated in the 
inset in Fig. \7}p with all the other control parameters 
and initial values being the same. 
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V. ATOMIC DYNAMICAL FRACTALS 

In this section, we treat the atom-photon interac- 
tion in a high-Q cavity as a chaotic scattering problem. 
Let us consider the scheme of scattering of atoms by 

(a) 



atoms 




FIG. 8: (a) Schematic diagram showing scattering of 
atoms at the standing wave and (b) sample atomic tra- 
jectories. 

the standing wave shown in Fig. Atoms, one by 
one, are placed at the point x = Q with different initial 
values of the momentum pa along the cavity axis. For 
simplicity, we suppose that they have no momentum 
in the other directions (ID-geometry). We compute 
the time the atoms need to reach one of the detec- 
tors placed at the cavity mirrors. The dependence of 
this exit time T on the initial atomic momentum po is 
studied under the other initial conditions and param- 
eters being the same. To avoid complications that 
are not essential to the main theme of this section, 
we consider the cavity with only two standing-wave 
lengths. 

A. Fock fractal 

In this subsection the field is supposed to be initially 
prepared in a Fock state. Before injecting into a cav- 
ity, atoms are suppose to be prepared in the superposi- 
tion state with ii„_i(0) = u„(0) = Vn-i(0) = fn(0) = 



0, z„_i(0) = —1/2, z„(0) = 1/2, i. c. in the state with 
zero population inversion z(0) = z„-i(0) + z„(0) = 0. 

At exact resonance {5 — 0) with u„_i(t) = 
Un{T) = 0, the optical potential U = {y/nun-i + 
y/n+ 1 Un) cosx — 5(zn-i + Zn) /2 is cqual to zero, and 
the analytical expression for the dependence in ques- 
tion can be easily found to be the following: 

r((5 = 0)-37r/2apoifpo>0, 

T{5 = 0) = 7r/2apo if Po < 0. ^ ' 

Atoms simply fly through the cavity in one direc- 
tion with their initial constant velocity and are reg- 
istered by one of the detectors. Out off resonance 
{5 7^ 0), the atomic motion has been numerically 
found in preceding section to be chaotic with posi- 
tive values of the maximal Lyapunov exponent in the 
following ranges of the values of the control parame- 
ters: the detuning | (5 | < 2 and the recoil frequency 
a ~ IQ-'^ H- 10-2. Fig. H shows the function T{pq) 
with the normalized detuning 5 = 0.4, the recoil fre- 
quency a — IQ-'^, and the initial number of cavity 
photons n = 10. The exit-time function demonstrates 
an intermittency of smooth curves and complicated 
structures that cannot be resolved in principle, no 
matter how large the magnification factor. Fig. 
shows magnification of the function for the small in- 
terval 64.1 ^ Po ^ 64.6. Further magnification in the 
range 64.2743 < po 64.2754 shown in Fig.Et reveals 
a beautiful self-similar structure. Some structures in 
Fig-Et that look like fractal are not, in fact, unresolv- 
able and self-similar. Magnification of the structure 
in the range 73.2 ^ po ^ 73.8 demonstrates quite a 
smooth function without unrcsolvable substructures 
and with only two singular points on the borders of 
the respective momentum interval. Beating in all the 
structures of the atomic fractal in Fig.^should be at- 
tributed to the structure of the Hamilton-Schrodinger 
equations l|5n|l which describe two atom-field oscilla- 
tors with slightly different frequencies. 

The exit time T, corresponding to both smooth and 
unresolved po intervals, increases in average with in- 
creasing the magnification factor. It follows that there 
exist atoms never reaching the detectors in spite of 
the fact that they have no obvious energy restrictions 
to leave the cavity. Tiny interplay between chaotic 
external and internal dynamics prevents these atoms 
from leaving the cavity. The similar phenomenon in 
Hamiltonian systems is known as dynamical trapping 
pgf . Different kinds of atomic trajectories, which 
are computed with the system ()30|l . are shown in 
Fig. ^p. A trajectory with the number m transverses 
the central node of the standing- wave, before being de- 
tected, m times and is called m-th trajectory. There 
are also special separatrix-like mS'-trajectories follow- 
ing which atoms in infinite time reach the stationary 
points Xs — ±7rs (s = 0, 1, 2, . . . ), ps = 0, transvers- 
ing m times the central node. These points are the 
anti-nodes of the standing wave where the force acting 
on atoms is zero. A detuned atom can asymptotically 
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FIG. 9: Fock atomic fractal with different resolutions. 

reach one of the stationary points after transversing 
the central node m times. The trajectory with number 
1, showing in Fig. 03, is close to a separatrix-like IS- 
trajectory. The smooth po intervals in the first-order 
structure in Fig. ^ correspond to atoms transversing 
once the central node and reaching the right detector. 
The unresolved singular points in the first-order struc- 
ture with T — oo at the border between the smooth 
and unresolved pq intervals are generated by the 15- 
trajectories. Analogously, the smooth po intervals in 



the second-order structure in Fig. ^jp correspond to 
the 2-nd order trajectories with singular points be- 
tween them corresponding to the 2S'-trajectories and 
so on. 

There are two different mechanisms of generation of 
infinite exit times, namely, dynamical trapping with 
infinite oscillations (m — oo) in a cavity and the 
separatrix-like motion (m ^ oo). The set of all ini- 
tial momenta generating the separatrix-like trajecto- 
ries is a countable fractal. Each point in the set can 
be specified as a vector in a Hilbert space with m in- 
teger nonzero components. One is able to prescribe to 
any unresolved interval of a m-th order structure a set 
with m integers, where the first integer is a number 
of a second-order structure to which trajectory under 
consideration belongs in the first-order structure, the 
second integer is a number of a third-order structure 
in the second-order structure mentioned above, and so 
on. Such a number set is analogous to a directory tree 
address: "<a subdirectory of the root directory>/<a 
subdirectory of the 2-nd level>/<a subdirectory of 
the 3-rd level>/...". Unlike the separatrix fractal, the 
set of all initial atomic momenta leading to dynam- 
ically trapped atoms with m ~ oo seems to be un- 
countable. 

Nonlinear fractal dynamics implies specific statisti- 
cal properties of chaotic motion in Hamiltonian sys- 
tems (for a review see \2^). Even a finite-dimensional 
atom-field system's phase space has a complicated 
topology whose simplified image is given by projec- 
tions of Poincare sections and fractal scattering func- 
tions. Tiny interplay between all the degrees of free- 
dom is responsible for dynamical trapping of atoms 
even in a very short microcavity and for anomalous 
statistical properties of the fundamental atom-photon 
interaction. The probability distribution of exit times 
P{T) for 2 • 10^ events with initial atomic momenta 
in the range 8 < po ^ 40 is shown in Fig. llOb in 
double logarithmic scale. It is close to a Poissonian 
distribution with comparatively short exit times up 
to T ~ 300 and demonstrates an algebraic decay 
with the characteristic exponent 7 ~ —3.72 at its tail 
(< 300T < 4000). In the range of initial atomic mo- 
menta 40 < Po < 41, the respective PDF for 1.5 • 10^ 
events is computed under the same other conditions 
to be a Poissonian-like. 



B. Coherent fractal 

Following to the scheme in Fig. let us consider 
in this section scattering of atoms by the standing 
wave field initially prepared in the coherent state H19|) . 
Before injecting into a cavity, atoms supposed to be 
prepared in the excited state |2) with 2;(0) = 1. In 
this case the initial conditions are 

Mn(0) = t^„(0)-0, Vn, 
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FIG. 10; Exit-time distributions in the Fock field: (a) an 
algebraic decay in the range of initial atomic momenta 
8 < po < 40 and an exponential decay in the range of 
initial atomic momenta 40 < po < 41. 



Fig. llll demonstrates the respective exit-time function 
T{pq) with S = 0.1, a = 0.001, and < n >= 10. A 
coherent state is not a sharply defined state, as a Fock 
one, but a superposition of an infinite number of Fock 
states. We used a truncated basis of 1000 Fock states 
for the cavity mode in our simulation. Resolution of 
one of the unresolved structures in Fig. Illb is shown 
in Fig. Illb . It resembles the respective fragment of 
the Fock fractal in Fig. ^jp. Further magnification of 
the function in Fig. Illb . which is shown in Fig. Illb . 
reveals again a self-similar structure with smooth and 
singular zones. 

We collect an exit time statistics with 6 • 10^ events 
by counting atoms with initial momenta in the range 
9 < Po < 30 reaching the detectors. The plot of the 
respective histogram of exit times, shown in Fig. 1121 
demonstrates a few local maxima. The corresponding 
times of exit of atoms can be approximately estimated 
with the help of the formula H35|) for resonant atoms. 
The first maximum around T ~ 150-^160 corresponds, 
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FIG. 11: Coherent atomic fractal with different resolu- 
tions. 



mainly, to atoms which, being initially directed to the 
right with comparatively low momenta pq ~ 10, turn 
back before reaching the central node and are regis- 
tered by the left detector. It follows from H35|l that 
for such atoms Ti ~ 7r/2apo — 157 at apo = 0.01. 
However, the atoms, to be injected initially with the 
momentum po — 20 25 and registered by the right 
detector, may contribute to the first maximum as well 
because, after transversing the central node, they can 
be accelerated and gain the values of the momentum p 
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FIG. 12: Exit-time distribution in the coherent field. 
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up to 50. The second local maximum around T ~ 200 
corresponds, mainly, to atoms with pq ~ 20-;- 30 which 
transverse the central node and are registered by the 
right detector for the time T2 ^ 3Tr/2apQ. The other 
local maxima of the PDF are not so pronounced as the 
first ones, they are formed by the atoms transversing 
the central node a few times. The PDF in Fig. [T^ 
demonstrates an exponential decay at the tail up to 
the exit times T 700. 



C. Bose-Einstein fractal 

In conclusion of this section we present in Fig. 13 
the fractal computed with the field initially prepared 
in a Bose-Einstein state (|20ll and atoms supposed to 
be prepared in the excited state |2) with z(0) = 1. 
The initial conditions are 



Un{0) = f„(0) = 0, 
Zn[0) = 



Vn, 



(36) 



(1+ < n >)"+!■ 



The values of the control parameters are the same as 
for the coherent fractal. Fig. 13 demonstrates again a 
self-similar structure of the exit-time distribution. We 
want to stress that the Hamilton-Schrodinger sets of 
equations generating the Fock and the other fractals 
differ strongly in the number of equations. Neverthe- 
less, the respective T(pq) functions are rather similar 
in their main features. The atomic fractals we have 
found are generated by dynamical chaos in the atom- 
photon interaction but not by noise. 



VI. CONCLUSION 

We have studied in the strong-coupling regime the 
fundamental interaction between a two-level atom 
with recoil and a quantized radiation field in a single- 
mode cavity in the mixed quantum-classical formalism 
modelling quantum evolution of the electronic-field 
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FIG. 13: Bose-Einstein atomic fractal with different reso- 
lutions. 



purely quantum system to be disturbed by transla- 
tional atomic motion. We have managed to derive 
an infinite set of the Hamilton-Schrodinger nonlinear 
ode's for real valued electronic-field probability am- 
plitudes and expectation values of the atomic center- 
of-mass position and momentum. These equations of 
motion are able to demonstrate known features of the 
atom-field quantum evolution in the strong-coupling 
limit such as collapses and revivals of the atomic popu- 
lation inversion modulated by atomic motion through 
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the nodes and antinodes of the standing wave. We 
have shown that even in the absence of any other inter- 
action with environment the Hamiltion-Schrodinger 
dynamical system provides the emergence of classical 
Hamiltonian dynamical chaos from cavity quantum 
electrodynamics . 

Wc have investigated the atom-photon nonlinear 
dynamics with Fock, coherent and Bose-Einstein 
quantum states of the initial cavity field. Positive val- 
\ies of the maximal Lyapunov exponent A have been 
found with reasonable values of the control parame- 
ters, the detuning of the atom-field resonance, 6, the 
atomic recoil frequency, a, and the initial mean num- 
ber of photons. Exponential sensitivity to initial con- 
ditions may manifest itself in the chaotic dependence 
of the output values of the atomic population inversion 
on its input values that may, in principle, be measured 
in real experiments. New manifestations of Hamilto- 
nian chaos in the atom-photon interaction, we have 
found numerically in the dependence of the atomic 
exit times on the initial momentum, are atomic frac- 
tals. They were classifield by the names of the respec- 
tive quantum field states, i.e., the Fock, coherent and 



Bose-Einstein fractals. We have found anomalous sta- 
tistical properties of the chaotic atom-photon interac- 
tion: Levy atomic flights through a cavity and a power 
law at the tail of the exit-time distribution function 
in the Fock field. Power-law tails should appear (at 
least, in some ranges of the control parameters) in the 
PDF's with coherent and Bose-Einstein fields, as well, 
but their detection would require a longer computa- 
tion time than we have used. 
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